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ABSTRACT 

The Wisconsin far-ultraviolet spectrometer aboard 
OAO-2 observed the wavelength region near 1216 A for 
69 stars of spectral type B2 or earlier. From the 
strength of the observed interstellar L a absorption, 
atomic hydrogen column densities were derived over 
distances averaging 300 pc away from the sun. Mea- 
surements of the equivalent widths of this absorp- 
tion were found to be inappropriate for deducing 
column densities because of the blending of nearby 
strong stellar lines by the 12 A wide instrumental 
profile of the OAO. Instead, the OAO data were com- 
pared to synthetic ultraviolet spectra, originally 
derived from earlier higher resolution rocket obser- 
vations, which were computer processed to simulate 
the effects of absorption by different amounts of 
hydrogen followed by the instrumental blending. 

An average volume density njj = 0.6 atoms cm -3 was 
derived using all the stars. For stars nearer the 
sun (r<140 pc) we obtained nn = 0.25 atoms cm -3 . 

The fact that no stars were observed much beyond 
1 kpc precludes our obtaining any general informa- 
tion on large-scale galactic structure. However 
there is evidence for a pronounced enrichment of 
gas toward a number of stars in Scorpius, and a 
general deficiency seems to exist in the 180° < l 11 
< 270° sector. A comparison of L a column densities 
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to those derived from observations of 21-cm emission 
shows a very poor correlation, which demonstrates 
the importance of the unavoidable differences in 
sampling geometry. The L a column densities corre- 
lated reasonably well with D line, H and K line and 
color excess measurements for the respective stars, 
and the average ratios <C N Nal/ N HI^ = x , 

j/Ng =2.5 x 10“ 9 and ^Ng/E(B— V)^ = 5 x 

10 21 atoms cm -2 mag -1 were obtained. For assumed 
ionization equilibrium constants and a model for 
the interstellar medium, we find the observed ratios 
imply the total abundances <Na/H> = 2.3 x 10“ 7 and 
<Ca/H> = 6.8 x 10 -9 . For a very recent grain model 
we derive <^ p Hl/ P dust> ~ 100 f rom the <N H /E(B-V)> 
relation. 


I. INTRODUCTION 

It has long been recognized that ultraviolet observations 
from space observatories can provide important new information 
on the interstellar gas and dust (Spitzer and Zabriskie 1959; 
Code 1960) . The successful launch and operation of the second 
Orbiting Astronomical Observatory (OAO-2) (Code, Houck, McNall, 
Bless and Lillie 1970) has demonstrated that stellar space ob- 
servatories can be designed to yield high-quality astronomical 
data for long periods of time. OAO-2 was launched on December 
7, 1968, and at the time of writing (~2-l/2 years later) con- 
tinues to collect high-quality scientific data. In this paper 
we will report on observations of the absorption from the 
interstellar Lyman alpha (L a ) line of atomic hydrogen at 1215.7 
A made by OAO-2. Other papers have dealt with OAO-2 observa- 
tions of interstellar dust (Bless and Savage 1970, 1972). 

Interstellar L a absorption was first observed with sounding- 
rocket instrumentation by Morton (1967) in the far ultraviolet 
spectra of 0 and B type stars in Orion (see also Jenkins and 
Morton 1967). Since that first observation, sounding rockets 
have recorded the L a line in approximately 20 stars. Reviews 
of these rocket observations and their implications have been 
given by Carruthers (1970a) and Jenkins (1970a). The primary 
conclusion which has been drawn from these observations is 
that the local interstellar hydrogen density is of the order 
of 0.1 atom cm -3 instead of the 0.7 atom cm - 3 generally quoted 
by 21-cm observers (Kerr and Westerhout 1965) . In the case of 
the rocket data the coverage of the sky has been small and 
therefore these early observations can not be considered as a 
general survey of the neutral hydrogen near the sun. However, 
these data have been very important because they have provided 
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an independent means of obtaining information on neutral inter- 
stellar hydrogen . 

The short wavelength spectrometer in the University of Wis- 
consin experimental package on board 0A0-2 has observed a 
large number of stars of type B2 or earlier in the spectral 
region 10’50-1800 A... For such stars the stellar L a line is in 
most cases narrow enough to reside within the highly saturated 
care of the interstellar L a profile. The distribution in ga- 
lactic latitude and 1 longitude of the stars observed is shown 
in Figure' Most stars are within 1000 pc of the Sun and are 
located near the plane of the galaxy. These new data provide 
a fairly representative sample of the local interstellar 
hydrogen ., 

A preliminary discussion of the OAO L a data was given by 
Savage and. Code (19701 ., In that analysis it was recognized 
that the relatively low resolution of the spectral scans ob- 
tained. by the OAO produced a blending of the interstellar L a 
line with, adjacent stellar lines,.. The equivalent widths re- 
ported by Savage and Code were for the blend of stellar and 
interstellar lines near 1216 A. This blend equivalent width 
was used, to obtain an upper limit to the hydrogen column den- 
sity-. For many stars the OAO blended equivalent widths were 
significantly larger (2 to 3 times) than the reported rocket 
La equivalent widths. In contrast, for other cases there was 
a fair agreement between the OAO blended equivalent widths and 
the rocket L a equivalent widths. Since that preliminary OAO 
report, a large effort has gone toward trying to understand 
the apparent differences between the OAO data and the rocket 
data and toward developing a method which corrects for the line 
blending in the OAO data. In this paper we present the re- 
sults of this investigation. 

II. OBSERVATIONS 

The OAO L a observations were obtained with an objective 
grating scanning spectrometer having an effective collecting 
area of 265 cm 2 . Starlight was incident on a plane reflect- 
ing grating with 300 lines/mm and the first order dispersed, 
light was focused by a parabolic mirror back through a hole in 
the center of the grating. Behind the grating, a slit having 
a width of 10 A was positioned in front of an EMR 541F' photo- 
multiplier which contained a LiF window. The spectrum was. 
scanned by a rotation of the grating in discrete wavelength 
steps of approximately 10 A. The integration time at each 
step position for all the L a data reported, here was 8 seconds. 
In this integration time an unreddened 3rd. magnitude BO star 
produces approximately 8000 counts in the region of 1200 A. 
Typically it was possible to obtain two spectral scans of a 
star from 1000-2000 A during one night period of an orbit. 
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All observations reported here were made at night since the 
presence of scattered sunlight made the daytime observational 
material difficult to interpret. For a more detailed descrip- 
tion of the OAO instrumentation and observational problems see 
Code et al. (1970) . 

Figure 1 shows some typical OAO-2 scans over the wavelengths 
1050-1400 A. We have plotted spectrometer counts vs. wave- 
length. An individual scan of an early type star would have 
the appearance of the plot for tt Sco (B1 V + B2, m v = 2.92). 

In this scan one can easily see spectral features that are pri- 
marily due to C III (1175.7 A, a multiplet of 6 lines), H I 
(1215.7 A), and Si IV (1393.8 and 1402.8 A). In addition, 
some weaker lines are apparent at 1300 A and 1335 A. The line 
at 1300 A is produced by a multiplet from an excited level of 
Si III (lines between 1294.5 and 1303.3 A) combined with the 
weak 0 I interstellar line at 1302.2 A. The line at 1335 A is 
due to a combination of stellar and interstellar C II (1334.5 
and 1335.7 A) . These line identifications were made with ref- 
erence to higher resolution rocket spectra of the Princeton 
group (Bohlin 1970) . 

The wavelength scale for the OAO data was established 
through a comparison of OAO data and higher resolution rocket 
data degraded to the OAO instrumental resolution. The pre- 
sence of the strong stellar lines makes such a comparison pos- 
sible (see Figure 1) . It was found that the number of Ang- 
stroms per spectrometer step depended on the direction of the 
scan (because of mechanical backlash) , the wavelength and the 
incident angle. 

When more than one spectral scan was obtained for an object, 
a detailed spectrum could be produced by combining scans. In 
Figure 1 we have combined two scans of the star x Sco (B0 V, 
m v = 2.83) and also four scans of the star <5 Sco (B0. 5 IV, m v 
= 2.32). These observations were obtained by utilizing the 
satellite's bore-sighted startracker. With this startracker 
it was possible to make a scan, change the incident angle by 
0.25 arc min (which corresponds approximately to a 2.5 A shift 
in the spectrum) , and then make another scan. If this process 
were continued one could produce a composite spectral scan 
similar to that for <5 Sco shown in Figure 1. However, only a 
few stars were observed using this technique to fill in the 
data between the 10 A stepping interval of the scanner. In 
general we combined two or more totally independent scans and 
shifted them in wavelength until a reasonable wavelength 
alignment was obtained. In fortunate situations this proce- 
dure led to a significant gain in information about the de- 
tailed shape of the spectrum. In other cases the data points 
simply all coincided. 

The OAO spectrometer records a background signal which is 
a combination of thermal dark counts, counts due to charged 
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F Igure 1. — 0A0-2 spectrometer scans o f the stars $ Sco 

(BO. 5 II/), t Sco (BO V) and tt Sco (87 V + B2}.« We have 
plotted n.aw digital counts vs. wavelength. The 1164 multi- 
plier. { \or the digital counts Is due to a scaler. In the 0A0 
electronics . The plot for tt Sco Illustrates an Individual 
scan. The plots lor x Sco and 6 Sco are a composite of 2 
and 4 scans , respectively. M ore frequent sampling than the 
spectrometer's 10 ft stepping Interval was accomplished by 
controlled onsets In Incident angle provided by the space- 
craft's bore-sighted startracker. The dashed line repre- 
sents the extrapolated background level. 
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energetic particles, sky background light, and scattered light 
in the optical system. The background due to charged parti- 
cles from the radiation belts is usually unimportant (for 8- 
second integration times) except at those times when the satel- 
lite passes through the South Atlantic radiation anomaly. In 
general, astronomical observations were not carried out during 
this period of time. The scattered light contribution domi- 
nates over the other sources of background for 0 and B type 
stars brighter than approximately m v = 2 . The scattering is 
undoubtedly from the grating and dust elsewhere in the optical 
system. The sky background light is mainly resonance scatter- 
ed L a from hydrogen in the earth's geocorona. The net back- 
ground due to all the above mentioned effects can be measured 
by observing an object shortward of the LiF cutoff of the pho- 
tomultiplier. In the data illustrated in Figure 1 one can see 
that for A.<1050 A (the LiF cutoff) the counts level off to a 
constant background level. For the purpose of obtaining L a 
line profiles we have assumed that the background level mea- 
sured at A<1050 A can be extrapolated as a constant background 
level for wavelengths greater than 1050 A (see dashed lines in 
Figure 1) . This assumption would have been invalid if some of 
the various background contributions were wavelength dependent 
or time dependent. A number of spectral scans of just the 
dark sky existed which demonstrate that the sky background is 
essentially constant with wavelength. We feel it is reason- 
able to assume that the scattered light contribution to the 
measured background is also wavelength independent. Since for 
most of the observational material presented in this paper the 
background is small compared to the observed signal at 1200 A, 
the relative uncertainty in the derived intensity will be neg- 
ligible. However, for stars faint in the ultraviolet the back- 
ground can be comparable in size to the signal being measured 
at 1200 A. Under this situation an uncertainty in the back- 
ground correction can produce a large uncertainty in the de- 
rived intensity. It is this uncertain background correction 
for faint objects and not counting statistics that limit the 
instrument for L a studies to 0 and B type stars brighter than 
about m v = 5. 

The data in Figure 1 illustrate the reproducibility of the 
OAO spectrometer observations. In general, spectral scans 
agree to within an uncertainty governed by photon statistics. 
Spectral scans of the same objects made two years apart have 
been compared, and typically they are identical to within 2 to 

O 9- 

o • 

Much of the analysis which follows relies upon a detailed 
fitting of the OAO scanner data to reconstructions of proto- 
type stellar spectra whose resolutions have been degraded. To 
accomplish this comparison one must know the shape of the 
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scanner's instrumental slit profile. This profile was derived 
from OAO observations of zero-order star images and from spec- 
tral scans of the relatively narrow He II (1640.4 A) emission 
line in the star Z Pup (05 f, m v = 2.25). Rocket observations 
of this line by one of us (EBJ) show its intrinsic width to be 
approximately 4 A. The zero-order image data were used to 
establish the shape of the instrumental profile while observa- 
tions of the He II emission line were used to determine width 
of the profile. A small correction was made for the finite 
width of the He II line. The derived function has an approxi- 
mate Gaussian shape with a full width at half intensity of 
12 A. When this function was used to degrade higher resolu- 
tion rocket spectra, a good agreement was obtained between the 
OAO scans and degraded rocket observations (see Figure 2) . 
These comparisons will be discussed in detail in § III. 

III. DERIVATION OF COLUMN DENSITIES 
a) Preliminary Considerations 

As mentioned earlier, some equivalent widths reported by- 
Savage and Code (1970) differed significantly from those of 
the earlier rocket investigators. It was not clear at that 
time to what degree errors in either case could be attributed 
to (a) faulty measurements, (b) an improper estimate of the 
stellar continuum level on either side of 1216 A or (c) con- 
fusion of the L a profile with nearby stellar absorptions. We 
shall briefly review our conclusions on these possibilities 
and show how they influence our approach in the reduction of 
the OAO data presented here. 

In spite of the measurement problems discussed in the pre- 
vious section, the photometric accuracy of the OAO photoelec- 
tric scans is generally more trustworthy than that of the 
rocket spectra recorded on film. Uncertainties from film 
grain noise, errors in the film's characteristic curve and 
ambiguous intensities below the film's threshold must all be 
recognized as possible sources of error in the raw data from 
many of the rocket observations. On the other hand, the rock- 
et spectra clearly resolve wavelength intervals significantly 
smaller than the width of the L a profile, and we depend upon 
our ability to see this finer structure to answer questions 
about the possible sources of error (b) and (c) mentioned 
earlier. 

A necessary (but not sufficient) test of the reliability of 
the OAO and rocket data is a comparison of the two after the 
rocket spectra have been convolved with the OAO instrumental 
profile. It is encouraging to see in Figure 2 that there is 
a surprisingly good agreement between the independently de- 
rived spectra. Figure 2 shows the references where the origi- 
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figure 2. — A compafii* on ofi 0A0 data (dot-5) with htgh n. nota- 
tion socket *pecVta which have, been smoothed by the 0A0 In 
6tn.urm.ntaZ pao^ile (cusive*). Both t\ te 0A0 and socket *pec 
tfia have been Keno fimalized to make the apparent continuum 
level nelatively {, lat at 1.0 (| \oa the socket *pectn.um o {, <$ 
Sco, thi-5 operation i-5 demon* tKated in the change &n.om Pig 
4a to 4b) . 
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nal rocket spectra were first recorded. The most serious dif- 
ference between rocket data and 0A0 data shown in this figure 
is seen in C Pup. In this case both the emission and absorp- 
tion components of N V look somewhat different, although the 
match in the core of the L a line is good. For this star the 
rocket intensity calibration may be somewhat in error for the 
strong emission of the N V line. The slight difference in the 
core of the La line for the y Ori comparison can probably be 
understood in terms of too much smoothing of the rocket obser- 
vation. In the convolution of the rocket data with the 0A0 
slit function we did not allow for the finite resolution of 
the rocket data. The rocket spectrum of y Ori was of slightly 
lower resolution (~3 A) than the other rocket spectra used. 

In the case of 6 Sco and tt Sco it appears that atmospheric 02 
absorption could produce the difference between the 0A0 and 
rocket spectra at 1240 A. One would expect a slightly weaker 
feature also to occur at 1206 A (Watanabe, Zelikoff and Inn 
1953) . The only arbitrary adjustment which was applied to 
force a correspondence for the profiles in Figure 2 was a mul- 
tiplication of each spectrum by a smoothly varying function 
which made the overall height of the continuum nearly level 
over the complete wavelength span. This renormalization com- 
pensated for any sloping trend in the stellar output and dif- 
ferences in each of the instrumental sensitivities versus 
wavelength. An illustration of this procedure is shown in the 
transformation of the rocket spectrum of 6 Sco between Figures 
4a and 4b. 

Figure 2 makes it very clear that for these stars previous 
differences in column densities from OAO and rocket data can- 
not be attributed to any fundamental disagreement in the re- 
corded information. In trying to explain why many of the rock- 
et investigations revealed significantly lower column densities 
than their own, Savage and Code (1970) suggested the possibi- 
lity that their predecessors had failed to account for the de- 
pression of the true continuum level by the broad damping wings 
in the L a profile, and the rocket equivalent widths were mea- 
sured under continua drawn systematically too low. Later, 
Jenkins (1971) reanalyzed the Princeton rocket L a data for 6, 
e and £ Ori by evaluating the goodness of fit of the observed 
profiles with theoretical absorption profiles as the hydrogen 
density and continuum height were both allowed to vary as free 
parameters . Jenkins pointed out that unless a careful reduc- 
tion procedure is used, rocket column densities could be sub- 
ject to errors. Nonetheless, the fact that Jenkins still 
found much lower column densities for the Orion stars than 
those of Savage and Code, together with preparations for the 
work presented here, made it clear that when interstellar col- 
umn densities were not very large, a principal contribution to 
the measured equivalent widths in the OAO absorptions recorded 
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near 1215 A came from neighboring stellar lines blended in 
with the L a absorption. 

b) Derivation of N w 

The aforementioned difficulties associated with 0A0 equiva- 
lent width determinations forced us to consider alternate cri- 
teria for establishing hydrogen column densities. The follow- 
ing approach was adopted for the analysis of the 0A0 data. 

The stars shown in Figures 2 and 3, whose spectra were observed 
at good resolution by rockets, were selected as prototypes for 
the analysis of stars of identical or nearly the same spectral 
type observed by the 0A0. 

The actual column densities N a (atoms cm -2 ) to these stars 
were either taken from Jenkins' (1971) analysis or evaluated 
using a somewhat simpler technique suggested in his article. 

A reconstruction of the appearance of a spectrum in the ab- 
sence of any absorption by atomic hydrogen was produced by 
multiplying the observed spectrum by exp (+t) , where the opti- 
cal depth t equaled 4.26 x 10 -20 atom -1 cm 2 A 2 N a (X-A 0 ) -2 
(Jenkins 1970a, 1971). Near the line core at Xo , where exces- 
sive magnification of small errors would occur, the flux was 
assumed to be a straight line joining the reconstructed spec- 
trum on either side. Figure 3 shows each prototype spectrum 
with and without interstellar absorption. One can see in this 
figure that the nearest strong stellar lines, Si III (1206.5 A) 
and N V (1240 A) , are often a major source of line blending 
for the L a profile and contribute to varying degrees in the 
different spectral types exhibited. Occasionally these lines 
show significant negative velocity shifts. For 6 Sco (BO. 5 IV), 
7T Sco (B1 V + B2) and y Ori (B2 III) , Si III is the major 
source of blending, while for ? Oph (09.5 Vnn) , 6 Ori (09.5 
II-III) , ? Ori (09.5 lb) and £ Ori (BO la), both Si III and 
N V contribute to the line blending. N V is a very strong 
component in £ Pup (05 f) . 

The next step consisted of generating a set of synthetic 
0A0 spectra for these stars for various trial column densities 
Nf A substitution of different values for Nt into the equa- 
tion for x and multiplication of each spectrum by exp ( — t ) 
produced high resolution versions of the spectra with various 
amounts of absorption, and these were subsequently smoothed by 
the OAO slit function. As a final preparation of the curves 
for matching with the actual data from other stars, these syn- 
thetic curves were renormalized so that they all reached the 
same height in the vicinity of 1160 A and 1270 A. Figure 4 
illustrates the various stages described for one of the proto- 
type stars, 6 Sco (BO. 5 IV). 

Figure 5 illustrates how the prototype spectra were used to 
obtain the adopted hydrogen column densities Nh to each OAO 
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PI guae 3. — Eight high aeAolutlon [aenoamallzed] ao cket Apectaa 
choAen { $oa paoceAAlng -into prototype. pao^lleA i oa compaal- 
Aon with the OkO data. The uppea cuave man L a aepaeA entA 
the original IntenAlty [lowea cuave) multiplied by exp 
[t(W<x)]» wheat the beAt estimate \ [oa the actual column den- 
sity No. to the Ataa Ia Ahown In team o ^ 10 zo atom a cm -2 . 
The numbea In paaentheA eA Ia the code aAAlgned fioa deAlgna- 
tlng In column 12 o £ Table 1 which Ataa woa ua ed ^ oa the 
deolvatlon ofi My. 
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TlguAe 4. — The AtepA In pAepaAlng a prototype. ApectAum ^Aom 
an IntenAlty tAaclng ofi a Aocket ApectAum. The ApectAum 
Ia noAmallzed to the line Ahown tn [a) to give a level 
ApectAum In (6). k^teA multiplication by exp [x(W a J], the 
AeconAt Auction In (c) o & the AtellaA ApectAum In the ab- 
a ence o & L a abA oAptlon Ia deAlved. ThlA pAo $lle Ia multi- 
plied by exp [-t(W^)] to Almulate the appeaAance ofa the 
ApectAum {jOA vaAlouA column denAltleA bl£. The cuAveA In 
[d) aAe the AeAultA &OA the dl^eAent valueA ( 1 0 2 0 atomA 
cm -2 ) a^teA a convolution with the InAtAumental pAofille 
Ahown In the coAneA o 6 [d] . The pAofilleA aAe Aeady ft oa 
compaAlA o n with the 0A0 data a^teA AenoAmallzatlo n to a 
AtandaAd height In ( e) by a AtAalght line. 
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Figure 5. — An Must station o f the derivation of hydrogen col- 
umn densities t $or eight 0A0 s tars. The dots represent 0A0 
data Ion the star lifted to the left o f each profile. The 
solid line* are prototype spectra for various total column 
densities [marked In units o f 70 20 atom s cm -2 ). The proto- 
type star {Indicated In parentheses next to the lowest 
curve of each set) Is the closest match In spectral classi- 
fication to the 0A0 star. Below the spectral type and V o f 
each 0A0 star Is shown the adopted Ny and error estimate 
(A, 8 or C j . 
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star. The general method of reduction was to choose the pro- 
totype star which matched most closely in spectral type and 
luminosity class the OAO star being analyzed. One then ob- • 
tained the hydrogen column density by choosing the trial col- 
umn density that gave the best fit to the OAO data. For ex- 
ample, in the case of the star Ori (BO. 5 IV-V) shown in 
Figure 5, we used the prototype star 6 Sco (BO. 5 IV) to obtain 
the hydrogen column density. As is evident, a trial column 
density of about 10 21 atoms cm -2 applied to 5 Sco would give 
a good match with the OAO observations of c{) 1 Ori. 

Since in our analysis we employed as prototype spectra only 
the eight stars shown in Figures 2 and 3, it was not possible 
always to use a prototype star which matched precisely an OAO 
star's spectral and luminosity class. For example, in the 
case of 10 Lac (09 V) shown in Figure 5 we used the prototype 
star £ Oph (09.5 V) to obtain a column density of 7 x 10 20 
atoms cm -2 . We matched c Sco (B1 III) with the prototype 
tt Sco (B1 V + B2) to derive Nh = 2 x 10 21 atoms cm -2 . We 
shall discuss the possible errors introduced into the final 
column densities through these unavoidable differences shortly. 

The OAO hydrogen column densities as obtained by the above 
mentioned procedure are contained in Table 1. For each of the 
stars observed we have listed the following: (1) HD number, 

(2) name, (3) galactic longitude and (4) latitude (System II), 
(5) spectral type, (6) V magnitude, (7) B-V color, (8) E(B-V) 
color excess, (9) distance to star (in pc), (10) OAO hydrogen 
column density, (11) error and position code (see § IIIc and 
§ IVa) , (12) prototype star used to obtain the hydrogen column 

density (see Fig. 3), (13) mean interstellar space density njj, 

(14) 21-cm hydrogen column density Nh (see § IVa) , (15) x (see 

§ IVi) , (16) simplified model density (see § IVa), (17) Na I 

and (18) Ca II column densities (see § IVb) , and (19) central 
depth (in percent) of the diffuse feature at 4430 A. Wherever 
possible, photometric data for northern stars were taken from 
Iriarte et_ al_. (1965) and for southern stars from Cousins and 
Stoy (1963) . Spectral types are those given by Lesh (1968) or 
by Hiltner, Garrison and Schild (1969) for northern and south- 
ern objects, respectively. The determination of E(B-V) utili- 
zed the intrinsic colors of Johnson (1963) . Distances were 
calculated using the photometric data in the table and abso- 
lute visual magnitudes from Blaauw (1963) along with the re- 
lation: 

m v - M v = 5 log r - 5 + 3E (B-V) . (1) 

Although there are a number of special cases where other 
methods could establish more accurate distances, we decided to 
retain the same system of assigning distances for all our 
stars in order not to upset the basis for studying regional 
variations of the hydrogen column density and its correlation 
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(19) 

Ac (%) 
A4430 

3.0(2) 

4.1(2) 

3.8(2) 

3.1(5) 

7.2(2) 

3.4(2) 

0.2(5) 

1.0(2) 

2.0(2) 

(18) 

Ca II 
Column 
Density 
10 12 cm- 2 

3.2(13) 

2.0(13) 

<0.43(12) 

2.9(13) 

2.7(13) 

(17) 

Na I 
Column 
Density 
10 12 cm- 2 

CM CM x-N CM CM 

r-t H ^.N /-sN ^«9N ,-n T-* f-4 ^ <t <N ^ 

w *w Mf ^ H >9 H «H <H >» w w H H ^ H H 

H w w w ^ '-'•'-'w w *3 -» w ww www 

CM CM CT\ C"» ITV 00 \D 'O H • sr fMCM SO CO CM O <J\ 

O OCMcScMHO^-4»>-M cn OOO CN CM -I (N N 

V A A e V V .<? 

vO 

iCajsuaa 

ispow 

O 0\ <r f-i \0 O oo r»1 cm co vycO 

M M 1 C3 n M3 >0 00 H >0 ^ r— CM vD I'm. CO oO 

in 

H 

X 

0.72 

0.24 

0.05 

0.22 

0.60 

0.35 

0.65 

0.85 

0.89 

0.56 

0.77 

0.22 

0.63 

0.84 

0.60 

0.58 

+• 

«a- 

21 cm 
Hydrogen 
Column 
Density 
10 2 “ cm -2 

3.3(1) 

21.(3, 1,0) 

4.5(4) 
9.9(3, 0.6) 

9.9(3, 3.5) 

14.(3, 1.5) 

6.2(1) 
8.6(6) 
12.(3, 5.4) 

12.(1) 

18.(7) 

12.(4) 

11. (3, 2.7) 

18.(3, 2.3) 

21.(3, 4.7) 
21.(8, 1.4) 

18.(7) 
10.(4) 
18.(3, 1.9) 

8.5(1) 

10.(7) 

12. (3, 2.1) 

24.(3, 0.9) 

23. (8, 0.2) 
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Table 1. References 


(1) Goldstein and MacDonald (1969) (10) 

(2) Duke (1951) (11) 

(3) Takakubo and Van Woerden (1966) (12) 

(4) Hobbs (1969, 1971) (13) 

(5) Stoeckly and Dressier (1964) (14) 

(6) Habing (1968) 

(7) Menon (1958) (15) 

(8) Linblad (1966) (16) 

(9) Wampler (1966) (17) 

(18) 


Muller (1958) 

Clark (1965) 

Merrill et. al. (1937) 
Spitzer et al.. (1950) 
Routly and Spitzer 
(1952) 

McGee et. al. (1966) 
Dunham (1941) 

Howard et al . (1963) 

Herbig (1968) 


‘Measured from data of Carruthers (1969) by Savage and Code 
(19 70) . 

+For those 21 cm surveys where the beam was not centered on 
the star, we have listed after the reference number the 
angular distance (in degrees) between the star position and 
radio position. 


with other observables. References for the data on Na I, Ca 
II, 21-cm column densities, and 4430 A central depths are 
listed at the end of the table. When more than one prototype 
star is listed in column 12 we have evaluated column densities 
by using both stars and then averaging the results. Such a 
procedure was necessary for certain intermediate spectral 
types . 


c) Error Estimate 

For each column density listed in Table 1 we have indicated 
the uncertainty by a letter classification system. Classifi- 
cation A implies the probable error is such that the true Njj 
may be 30% lower or 40% higher than the quoted value. Class B 
corresponds to -40% and +60%, while C is -50% and +100%. For 
example, the star <t> X Ori (see Fig. 5 and Table 1) has an esti- 
mated error of class A and therefore has a column density of 
10 x 10 20 + 4 x 10 20 - 3 x 10 20 atoms cm -2 . These error esti- 
mates are based on a consideration of the quality of each OAO 
spectrometer scan and the degree of confidence we have in the 
use of the prototype spectra to reduce a particular OAO star. 

The varying degrees of uncertainty in the OAO column densi- 
ties arise from several sources. In a number of cases the OAO 
line profiles are poorly determined because of a large and un- 
certain background correction. In these cases an estimate of 
possible errors has been made by allowing the background cor- 
rection to vary over permissible limits and determining the 
corresponding hydrogen column density for each background se- 
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lected. A source of error common to all the OAO column densi- 
ties is that of an incorrect allowance for line blending. 

Errors of this type can occur for a number of reasons. For 
example, the stellar features in the prototype stars may dif- 
fer somewhat from the features in the OAO star being analyzed 
because of (1) a spectral type mismatch, (2) a luminosity 
class mismatch, (3) large differences in element abundances 
between the prototype star and OAO star, (4) other differences 
due to the peculiar nature of certain stars being studied, and 
(5) possible atmospheric O 2 absorption at 1206 A in some of 
the prototype spectra. We estimate the O 2 feature to have an 
equivalent width no greater than 1 A and hence this possible 
feature would be relatively unimportant compared with the 
much stronger Si III stellar absorption at the same wavelength. 
We have 'obtained information on how large errors of type (1) 
and (2) can be by using several prototype stars to analyze a 
single OAO star. Our general finding was that when the inter- 
stellar hydrogen column densities are large (Njj> 2 x 10 20 atoms 
cm -2 ) the errors due to spectral class or luminosity class 
mismatches are small (class A and B errors) . This was the 
consequence of line blending being less important when the 
interstellar L a line dominates the stellar lines. In contrast, 
when the interstellar hydrogen column densities are small 
(Nh< 2 x 10 20 atoms cm -2 ) the effects of line blending are more 
important and small mismatches . between the OAO star and the 
prototype star produce larger errors (generally class B and C) . 
For this reason the stars with small hydrogen column densities 
in Table 1 are usually of class C error while the stars with 
large column densities are usually of class A error except in 
those cases where the OAO profiles are poorly determined. 

As we saw in Figure 1 , for every star observed at L a the 
OAO spectrum scanner obtains information on a number of strong 
ultraviolet stellar lines, and we can verify that the prototype 
star with the most similar ultraviolet spectrum has been selec- 
ted. Since the Si III (1206.5 A) feature is a major source of 
line blending with L a in the OAO data, it is particularly in- 
structive to study the relative strengths of the prominent Si 

III (1300 A) and Si IV (1400 A) multiplets . Differences in 

these two silicon features between the prototype and subject 
star provides a gauge of the importance of any of the four 
earlier mentioned dissimilarities in reducing the accuracy of 
our L a absorption measurement. For example, in Figure 1 we 
can see that the stellar lines of Si III and Si IV in x Sco 

(B0 V) and 6 Sco (B0. 5 IV) are of similar strength. This ob- 

servation gives one confidence that the use of the rocket data 
on 8 Sco as a prototype for the column density determination 
in x Sco is valid. Through comparisons of this type we have 
noted that for stars of spectral type B0 to B2, the spectra of 
stars of luminosity class V, IV, III and II are quite similar. 
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while the supergiants have much stronger ultraviolet lines 
(Bless and Code 1970) . We therefore feel we are justified in 
using, for example, the prototype star y Ori (B2 III) to ana- 
lyze other B2 stars of luminosity classes ranging from V 
through II. For the O type stars the differences among lumi- 
nosity classes are larger. However, in the case of 09.5 stars 
we have prototypes for luminosity classes V, II and I. 

The column densities reported in Table 1 should in most 
cases pertain to the interstellar medium. Jenkins (1970a) has 
considered the possibility that L a emission from the star or 
its surroundings might significantly modify derived interstel- 
lar column densities. His conclusion was that such modifica- 
tions are very unlikely. A more relevant problem for the data 
presented here is that of determining if photospheric absorp- 
tion at L a is contributing to the observed hydrogen column 
densities. As we indicated earlier, stars of spectral types 
later than B2 were excluded from this paper because in such 
stars the stellar L a line is very strong. For all of the 
stars listed in Table 1, we have plotted Nh vs. spectral type 
and have found no overall tendency to obtain larger column 
densities in the cooler stars. It will be apparent from the 
material presented in § IVc on the gas to dust correlation .. 
that for B2 stars the stellar L a line does not significantly 
modify the derived interstellar column densities. On the 
other hand, a reasonable fraction of the L a absorption could 
be stellar when the measured Ny is as low as approximately 
10 20 atoms cm -2 ; hence it is probably wiser to consider these 
lower values as upper limits for the interstellar N H . 

d) Comparison of 0A0 and Rocket Column Densities 

A comparison of our results with various rocket column den- 
sities is given in Table 2. For the prototype stars we have 
listed the results obtained from measurements different from 
the ones we have employed for our comparisons. It is reassur- 
ing to note that the N a values we have used do not differ 
appreciably from the values of the other investigators. In 
the case of the prototype star y Ori, our determination of N a 
from Carruthers 1 (1968) data is somewhat larger than his own 

measurement. The Njj values for the remaining star£ in the 
right hand portion of the table were all derived from the OAO 
data by the technique discussed in § Illb. For these stars 
roughly half show good agreement (y Cas, x Ori, a Vir, B 1 Sco, 
a Sco and co 1 Sco) while the rest disagree considerably with 
the OAO column densities. For n Ori and a Ori the differences 
can be understood as due to the low quality of these particu- 
lar rocket spectra. For 8 CMa the OAO value has a C error 
estimate, and in addition, Carruthers' spectrum for this star 
was of low quality (private communication) . With k Ori, e Per 
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Table 2. Comparison of OAO and Rocket Column Densities 


Prototype Stars 


Star 

Spectral 

Type 

N a 

10 cm 

N H (rocket) 

Reference 

Y Ori 

B2 III 

0.65 

0. 39 + 5! • 19 
-0.1 5 

3 

6 Ori 

09.5 II-III 

1.3 

1.5 

1 

z, Ori 

09.5 lb 

1.75 

1.5 

1 




1 i +0.30 
' -0.26 

3 




1.1 +0 ' S . 
-0.4 

5 

S Pup 

05 f 

0.70 

0.59 +0 - H 
-0.16 

3 

IT SCO 

B1 V, B2 

6.3 

7 5 + ^ • 
/,b -4. 

6 

6 Sco 

B0. 5 IV 

12.5 

12 +6 • 
' -3. 

6 

? Oph 

09 . 5 Vnn 

5.0 

4.2 

10 




4 . 6± / . 5 
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3. Carruthers (1968) 7. Jenkins (1970b) 

4. Smith (1969) 8. Carruthers (1970b) 

5. Carruthers (1969) 9. Bohlin (1970) 

10. Smith and Stecher (1971) 
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Table 2, continued 
0A0 Stars 


Star 

Spectral 

Type 

N h (OAO) 
10 2 0 cm 2 

OAO 

Error 

Estimate 

Nh ( rocket) 
10 2 0 cm 2 

Reference 

Y Cas 

BO. 5 IVe 

1.9 

B 

1 . 0 +7 • 7 
-0.3 

9 

n Ori 

BO . 5 Vnn 

6.0 

A 

2.7 +3 - 3 

-2.0 

2 

0 Ori 

06p, 09 . 5 Vp 

11 . 

B 

4 8 +2 ‘ 7 

5 

a Ori 

09.5 V 

5.0 

A 

2 7 + 7 • 3 

Zm -1 .2 

2 

i Ori 

09 III 

1.5 

B 

1 5 + 7 * ^ 
X -0. % 

2 





n 7 * ® * 3 

3 





1 4 + 0.6 
* -0.4 

5 

k Ori 

BO. 5 la 

4.0 

A 

0 rq + 0.33 
°- 89 -0.2S 

3 

3 CMa 

B1 II-III 

0.2 

C 

0 57 + 9 * ^ 77 
U * -0. 34 

3 

a Vir 

B1 IV 

0.4 

B 

<0.4 

4 

3 1 Sco 

BO. 5 V 

13. 

© 

A 

14. ±6. 

7 

a Sco 

B1 III 

20 . 

A 

20 . ±3. 

7 

w 1 Sco 

B1 V 

12 . 

A 

16. ±4. 

7 

e Per 

BO. 5 III 

3.0 

B 

1 . 1± 0 . 3 

8 

£ Per 

07.5 

20 . 

A 

4.2 +7 - 9 
-/. 5 

8 
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and £ Per the column densities of Carruthers are all signifi- 
cantly lower than the OAO values . We believe the problem for 
these stars (and also his determination for Y Ori) could be 
related to that previously discussed by Savage and Code (1970) 
for 0 1 Ori C, for which it was shown that Carruthers (1969) 
had underestimated the equivalent width by placing a continuum 
level too low and by not allowing for the broad damping wings 
of L a . 


IV. DISCUSSION 

a) The Distribution of Hydrogen 

Of the 69 stars listed in Table 1, 81% are less than 500 pc 
away, and all but one are closer than 1 kpc. The proximity of 
observed stars is simply a manifestation of the magnitude lim- 
it for OAO observations of good quality at L a . Over such a 
small distance scale, the observed distribution of interstel- 
lar atomic hydrogen provides relatively little insight on gen- 
eral features in the structure of our galaxy, such as spiral 
arms. Moreover, the well known grouping of O and B type stars 
into associations makes it inevitable that some limited sec- 
tions of the sky are well sampled while other more extensive 
regions are left unobserved. In spite of these shortcomings, 
however, interesting conclusions on significant local excesses 
or deficiencies of gas from one region to the next may be 
drawn from the stars in our list. A strong suggestion of this 
variability came from earlier rocket observations (Jenkins, 
Morton and Matilsky 1969) . We have here a large enough set of 
observations, all made by the same instrument, to over-rule 
any earlier misgivings one may have had that the differences 
in hydrogen densities were attributable to systematic errors 
arising from the markedly different wavelength resolutions of 
the separate rocket flights or the differences in spectral 
types of the stars available. 

The information on distance r and average volume density n H 
in Table 1 is plotted in Figure 6a in the form of ellipses 
whose centers are located at the galactic coordinates for the 
appropriate star. The length of an ellipse's horizontal axis 
is proportional to r and the vertical axis corresponds to nfj. 

It follows that the area of each ellipse is proportional to 
the measured column density Njj . One can immediately see from 
this plot that the hydrogen gas distribution is not at all uni- 
form. Certain areas of the sky stand out as being characteris- 
tically rich or deficient in hydrogen. 

In trying to assess the degree of variability in the L a col- 
umn density measures, it is helpful to compare the observations 
with what one would expect to obtain if the hydrogen were per- 
fectly uniform except for a stratification parallel to the ga- 
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lactic plane. For the z dependence of the volume density njj 
we shall adopt the distribution given by Schmidt (1957) scaled 
to a thickness between half-density points of 220 pc (Schmidt 
1957, McGee and Milton 1964) . The density is normalized to 
0.71 atom cm -3 for z = 0, which corresponds to a compromise be- 
tween the Nh = 3.5 x 10 20 esc |b| atom cm -2 relation observed 
somewhat away from the galactic equator by McGee and Murray 
(1961) and the 2.5 x 10 20 esc |b| atom cm -2 result of Goldstein 
and MacDonald (1969) . The expected amount of hydrogen the 
model would give to each star is shown in column 16 of Table 1. 
The vertical axis of each ellipse in Figure 6 b has a length 
proportional to the L a N H , whereas the horizontal axis corre- 
sponds to the expected density for the simple distribution. 

The shape of the ellipses depict excesses or deficiencies of 
gas, which is basically the same information given in Figures 
4 and 5 of the article by McGee and Murray (1961) . 

With the exception of a few of the Orion stars, the observa- 
tions indicate a general deficiency of gas in the 180-270° sec- 
tor of galactic longitude (an average of nn over all of the 
stars in this area is 0.27 atom cm -3 ) . Three stars on the 
order of 1 kpc away and near the galactic plane, 139 Tau, 15 
Mon and UW CMa, have nn equal to 0.24, 0.061 and 0.12 atom 
cm -3 , respectively. This suggests that the lack of hydrogen 
in this area of the sky may not be confined to only the local 
few hundred pc distance which is typical of most of the stars 
covered in our survey. A deficiency in 21-cm emission near 
H = 240° may be seen in both the plots of integrated bright- 
ness temperature of McGee and Murray (1961) and the contours 
of peak brightness temperature shown by McGee, Murray and 
Milton (1963) (perhaps best seen in Figure 4 of the latter arti- 
cle) . This 21-cm minimum does not appear to extend as far back 
in longitude as 15 Mon ( Z = 203°) or 139 Tau (it = 184°), how- 
ever. 

Some of the reduction in Nh for 5 Pup (Z = 256°, b = -5°) 
may be attributed to the ionization of hydrogen within the Gum 
Nebula which surrounds this star (Gum 1955, Rodgers et al . 

1960) . The only other line of sight we observed which inter- 
sects this nebula is that toward UW CMa (Z = 238°, b = -5°), 
but this star is very distant, and only the outer edge of the 
nebula is in the way, so the ionization would probably be of 
negligible importance. Brandt, Stecher, Crawford and Maran 
(1971) have suggested that ionization of the interstellar gas 
occurs well outside the apparent boundary of the nebula. They 
made use of the earlier measured low Nh to £ Pup and y Vel, 
together with an estimate of 0.4 atom cm -3 for nn in the neigh- 
borhood of the sun, to argue that the volume of the ionization 
extended to within 60 pc of the sun. Such an interpretation 
would appear to be substantiated by our observation that neut- 
ral hydrogen seems to be deficient over a very wide angle in 
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this region. On the other hand, it is hard to understand why 
the dispersion measures of MP 0818 (S- = 237°, b = +12°), PSR 
0628-28 U = 237°, b = -17°) and CP 0834 U = 220°, b = +26°) 
are significantly less than PSR 0833-45, MP 0835 and MP 0736 
located well inside the nebula (Gott and Ostriker 1971) . It 
therefore seems more plausible that the ionization is still 
confined to the 40° diameter nebula, and that there is simply 
less interstellar gas in this particular area of the sky. 

Another area which is prominent in Figure 6a is a small 
patch in Scorpius where the gas density is particularly rich. 
Although there is a spur of 21-cm emission projecting out of 
the galactic plane in this general neighborhood, details of 
the emission contours mapped in Figure 1(a) by McGee et al . 
(1963) do not correlate closely with the variations from one 
Scorpius star to the next. This might be an indication that 
these stars are embedded in a cloud complex, with some stars 
in front and others behind or within the gas cloud. (The 
properties of the medium in this region have been discussed 
by Herbig (1968) and Johnson (1970, 1971)). 

It would be tempting to use the La results to actually de- 
termine the general behavior of njj(z) in the solar vicinity 
since in each case an observation of Nh sin |b| is an average 

; z (star) 

n H (z) dz along the line of sight. Unfortunately 

0 

the plot of log N H sin |b| vs. log z for our observations in 
Figure 7a shows considerable scatter, and it is difficult for 
us to define any systematic departures from the concept that 
n„ does not vary with z. (In this and succeeding figures in 
this paper, we have plotted the assigned A, B or C error code 
for each point instead of error bars to eliminate clutter.) 
Figure 7a, along with displays in Figure 6, helps to emphasize 
the strong irregularity in the H I gas distribution. We feel, 
however, that this type of display can be misleading because 
at high z we cannot avoid preferentially viewing toward isola- 
ted enrichments of gas projecting out of the galaxy, since it 
is precisely in such areas one usually finds young stars. If 
our sampling of directions were truly random (which it is not) , 
the plot would be more meaningful. Similarly, we feel that 
any attempt to find some simple distribution which best fitted 
our data (such as the determination of the plane of the local 
gas layer in the manner of Goldstein and MacDonald (1969)) 
will be more a reflection of the irregularity of our sampling 
procedure, rather than a true generalization on the local 
distribution of hydrogen gas. 

Figure 7a shows that 78% of the stars we observed have an 
average volume density nn somewhere between 0.1 and 1.0 atoms 
cm -3 along the lines of sight. An average nn for all of the 
stars equal to 0.6 atoms cm -3 was obtained over a mean dis- 
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tance of 300 pc from the sun. The average color excess per 
unit distance £ E(B-V)/E r for these stars is equal to 0.31 
mag kpc -1 . On the other hand, an average figure for the gen- 
eral reddening within 1 kpc of the sun is 0.61 mag kpc -1 
(Spitzer 1968) , which is roughly twice our average for the OAO 
stars. This difference demonstrates that our magnitude limit 
V£5 for L a scans imposes a bias against the more heavily ob- 
scured stars at large distances. Hence, if one were to assume 
that the ratio of hydrogen to dust is constant everywhere (see 
§ IVc and Fig. 9) , a representative figure for the density 
over all directions would be twice the average njj we have de- 
rived. There are seven stars closer than 140 pc for which the 
average nn is 0.25 atoms cm -3 . We do not expect these nearby 
stars to be subject to the aforementioned observational selec- 
tion. 

Figure 6a or 6b reveals that of all the stars observed, a 
considerable fraction share some characteristic with their 
neighbors. We have already singled out for discussion the hy- 
drogen-rich Scorpius region and a collection of stars having 
anomalously low column densities. The stars in Orion are all 
located within a reasonably confined solid angle, and stars in 
Crux, Centaurus and Lupus ( O ') , although fairly spread out, all 
have nearly the same distance modulus. It is of benefit to 
later discussions on the correlations of Njj with other observ- 
ables to define memberships in groupings, even though the 
assignments are somewhat arbitrary. We have designated four 
groups, shown by the dashed perimeters in Figure 6a, and have 
put the remaining stars into a classification of "other" since 
they seem to be reasonably heterogeneous with respect to Ny, r 
and position in the sky. The symbols accompanying the regions 
in Figure 6a identify the membership of each star in the corre- 
lation plots (Figures 7, 8 and 9) and column 11 of Table 1. 

In the introduction (§ I) we mentioned that the findings of 
the sounding rocket observers generally indicated the local 
hydrogen density was lower than previously suggested by over- 
all surveys of 21-cm emission. More specifically, flights by 
several investigators clearly demonstrated that only occasion- 
ally did one find a L a measurement yielding a column density 
as high as that suggested by the 21-cm emission in the same 
direction. This would not have been surprising for low galac- 
tic latitude stars since much of the 21-cm emission*would have 
originated from portions of the galaxy beyond the stars, but a 
number of stars with |b| £ 10° exhibited densities a factor of 
10 lower than the corresponding 21-cm measures (Jenkins 1970a) . 

We are now in a position to repeat the comparison with the 
21-cm data on a much wider scale. Figure 7b is a logarithmic 
plot of the L a Njj against the 21-cm measures listed in column 
14 of Table 1 for stars with |b| > 10°. For some of the stars 
we have observed, we could draw upon the results reported by 
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1.0 10 . 100 . 1000 
LOG Z(pc) — 


Figure 7 a . — % sin |fa| values plotted are an average measure a- 
long the llne-o fi-sxght o fa the. amount o ^ hydrogen situated be- 
tween thz z distance o a star and thz plane o & the galaxy 

/ z(s tar ) 

ny[z)dz li ny depends only upon z) . 

0 

Thz dashed line shows thz expected behavior o £ thz measure- 
ments ll ny[0) = 0,71 atoms cm ~ 3 with a decrease In density 
away firom the galactic plane allowing thz distribution 
given by Schmidt [1 957). 

Habing (1968), Goldstein and MacDonald (1969) or Hobbs (1971), 
where the 21-cm emission has been sampled from a small diame- 
ter (1/2°) beam centered on the star's position (for the ex- 
press purpose of comparing 21-cm and visual absorption line 
data). For many other cases, however, we had to rely on a 
measurement taken on the nearest grid point in the routine sur- 
veys of the other investigators referenced in the table. In 
each case we have shown in the tabulation the distance of the 
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F Igunelb . — A con nelatlon o i L a Ny measunements [at b> 10°) 

against 1 1 - cm emission meas uA.eme.nts [see column 14 In Tabic 
1) in the same direction. The lino, away in om the leit oi 
each symbol Indicate an estimated downwand connection oi 
the 21 -cm nesults to compensate ion backgno and emission be- 
yond the stans (see § I l/a] . The shapes oi the plotted 
points show the stan locations as defined In Fig. 6a, and 
the bolden symbols Indicate the small-beam obs envatlo ns 
dlnectly centened on the stans . 


survey point from the star. We expect to have the closest cor- 
respondence with the centered observations; for Figure 7b we 
have chosen these measurements, when available, in preference 
to the survey values and plotted the symbols more boldly to em- 
phasize their greater relevance to our measurements. For the 
surveys by Takakubo and van Woerden (1966) and Linblad (1966) 
we added up the N^ values given for the Gaussian components of 
optical depth, which were computed under the assumption that 
the spin temperature equalled 125°K. For all of the other ref- 
erences, we listed Ny simply computed from the integral of 
brightness temperature over frequency, which invokes the as- 
sumption that T 2l-cm << '*‘ * T ^ e f act that optical depths often 
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are not negligible means that the column densities computed by 
the latter method are usually somewhat lower than the true 
values. 

Figure 7b demonstrates that the OAO data show virtually no 
correlation with the 21-cm emission measurement. A large ma- 
jority of the column densities deduced from the 21-cm data are 
higher than those obtained from the L a absorption measurements. 
To some degree, one would generally expect to find higher Nh 
values from the 21-cm data because of the additional contribu- 
tion from hydrogen beyond the stars we have observed. To ob- 
tain a rough estimate of how much the 21-cm Nfj should be re- 
vised downwards to correct for the background emission, we may 
resort to the idealized model for nj^(z) discussed earlier and 
say it is reasonable to multiply each 21-cm result by the fac- 
tor 


x 


/• z (star) 
J 0 n H ( z ) dz 




0 n H (z)dz 


( 2 ) 


where the integral to infinity equals 3.0 x 10 20 esc |b| atoms 
cm -2 . The factor x is tabulated for each star in column 15 of 
Table 1, and each point plotted in Figure 7b has a line extend- 
ing toward the left (ending at an "x") whose length corresponds 
to the logarithm of the correction factor. Even with the x 
corrections, however, most of the points still lie well below 
the diagonal N H (L a ) = N^ (21-cm) line. For stars above this 
line, it is likely that saturation of the radio emission is 
causing us to underestimate the 21-cm N^. 


b) Correlation of H I with Na I and Ca II 


Measures of the interstellar D line absorption of Na I and 
the H and K lines of Ca II have been made for many of the 
stars in our program, as is evident from the many entries in 
columns 17 and 18 of Table 1. A comparison of H I to Na I and 
Ca II column densities should allow us to derive reasonably se- 
cure values for the overall abundance ratios in the inter- 
stellar medium, as well as to carry out a limited study on how 
much these ratios fluctuate over the various lines of sight. 
Past studies of this nature (Lawrence 1956, Howard, Wentzel 
and McGee 1963, Takakubo 1967, Habing 1968, Habing 1969, Gold- 
stein and MacDonald 1969, Hobbs 1971) have had to contend with 
the uncertainty of how well the optical column corresponded 
with the 21-cm beam. We have seen from the results of the 
previous section that even for small radio beams exactly cen- 
tered on a star, the 21-cm measurement is not very representa- 
tive of the column density to the star. The earlier investi- 
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gators, however, did have the opportunity to use the velocity 
information in both the 21-cm emission and visible absorption 
data to help resolve this problem. Interesting conclusions 
from the earlier investigations on the correlations of specif- 
ic velocity components cannot, of course, be verified or ex- 
plored further with the L a results. 

The equivalent widths of both components of the doublets 
have been obtained by the investigators referenced in the ap- 
propriate columns of Table 1, with the exception of Hobbs 
(1969) , who observed with an interferometric, photoelectric 
scanner and displayed high resolution (0.5 km sec -1 ) tracings 
of the D2 line (and also the Dl line when D2 appeared to be 
saturated) . For all but Hobbs ' data, we have made use of 
Stromgren's (1948) Table 2 to derive column densities by the 
doublet ratio method, except for a revision of the transition 
f values to those given by Wiese, Smith and Miles (1969) of 
0.327 for Dl and 0.344 for H. For the cases where the doublet 
ratio exceeded 2.0, we used the procedure of Routly and Spitz- 
er (1952) of equating W^i to (Wpl + Wd 2)/3 and assuming the 
lines were completely unsaturated. Occasionally, only an up- 
per limit for the equivalent width of one component was avail- 
able for a particular star. In this situation we derived col- 
umn density upper limits assuming the line was weak enough 
that no saturation occurred. The two entries for Ca II in 
6 1 Sco (HD 144217) exemplify that this assumption may not al- 
ways be valid and that such upper limits are not entirely 
trustworthy . 

The greater profile accuracy and detail in the sodium data 
of Hobbs (1969) warrant a treatment which does not rely upon 
the simplifying assumption on the velocity behavior which is 
inherent with the doublet ratio method. As long as the resi- 
dual intensity at the line center is not near zero, we can de- 
rive the neutral sodium column density from a straightforward 
integration of the measured optical depth over frequency, 

- i 

‘(km sec -1 ) cm 2 St dv (3) 


Optical depths t greater than 3.0 were hard to measure accur- 
ately and 4 . 6 was the maximum value we could include in the 
reduction. For stars having an absorption much stronger than 
x = 3 over a significant fraction of the profile's total ve- 
locity spread, the measurements very likely underestimate the 
true column densities, and hence these values were quoted as 
lower limits in Table 1. When available, Dl line values were 
listed in preference to the more saturated D2 line determina- 
tions. For some stars, we noticed that a large fraction of 
the t D 2 measurements ranged from 10 to 40% lower than the ex- 


N. 


me 


Na I -rre 2 f x 


/t dv = 


1 . 
9 . 


96 x 10 1 1 (Dl) 
8 x 10 1 0 (D2) 
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pected value of 2 t D i at the same velocity. This effect occur- 
red for T £)2 significantly less than 3, and we conclude this 
may be an indication that for some cases the 0.5 km sec -1 
resolution was not quite adequate for resolving saturated fine 
structure in the profiles. Hence it is possible that even 
some of the determinations not having excessive x values may 
slightly underestimate the true Na I column density. 

Occasionally, telluric water vapor absorptions appeared 
along with the D2 profiles. Hobbs pointed out this contamina- 
tion when it occurred in his data, and in such situations we 
either totally relied upon the Dl profile (if available) or we 
eliminated the water peaks from the integration if they were 
small and separable. At times we could make use of only D2 
tracings which had strong water lines blended with the sodium 
absorption, and in such cases we quoted the Na I column densi- 
ty values as upper limits in Table 1 and Figure 8a without 
attempting to separate out the spurious absorption since we 
could not estimate its strength. Stars having both large sodi- 
um optical depths and water vapor lines have values designated 
as approximate (M in the table. 

When the Na I column density has reached a value well in ex- 
cess of that necessary to saturate the D lines, the 3302 doub- 
let becomes strong enough to measure. Once again the ratio 
method can be used to derive a figure for N^ a j by now compar- 
ing the 3302 and D line strengths (Stromgren 1948) . Herbig 
(1968) used this approach to arrive at N Na j = 4.9 x 10 1 3 
atoms cm -2 to L, Oph (HD 149 757) . As noted by Herbig (196 8) 
and Hobbs (1971) , the doublet ratio method applied to the D 
lines fails for such large column densities, and it appears 
that even our value of 3.7 x 10 1 2 atoms cm -2 derived from in- 
tegration of the high resolution data is far lower than Herbies 
result. The discrepancy with our method is rather puzzling in 
view of the fact that Herbig 's curve of growth shows the opti- 
cal depth To = 56 for the center of the Dl line, whereas the 
largest x seen in the Dl tracing of Hobbs (1969) is about 2.4. 
The residual intensity measured at the bottom of the strong Dl 
profile seems, by a good margin, to be larger than either the 
statistical noise or any reasonable uncertainties in subtract- 
ing off the parasitic light and dark current contributions de- 
scribed by Hobbs, although the profile does appear to have a 
slightly flattened bottom which suggests saturation. 

In view of the aforementioned disagreement, we are somewhat 
skeptical of the larger N^ a i values quoted in Table 1. In 
deriving a representative average value for the density ratio 
<Na I/H I> we disfavored the larger NNa I measurements as 
well as those where only upper limits were available. Esti- 
mates of 3.5 x 10 -9 for <Na I/H I> and 2.5 x 10“ 9 for 
<Ca II/H I> were derived from a visual inspection of Figures 
8a and 8b, where points in the middle of the sequences were 
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given the greatest weight. Our value for <Na I/H I> is in 
excellent agreement, with an average’, given by Hobbs (19 71) of 
4 x 10- 9 . 

A large fraction of the interstellar sodium and calcium 
atoms are expected- to- be in the ionization state just higher 
than that, observed., since both Na I and Ca II can be ionized 
by starlight, photons less energetic than the Lyman limit at 
13.6 eV. Hence.,, to deduce representative values for the total 
element abundance ratios; <Na/H> and <Ca/H> , one. must eval- 
uate the equations for ionization equilibrium and accordingly 
apply large- corrections' to the observations of Na I and Ca II. 
In. the case of sodium, for any point along the line of sight 
we may consider the balance of photoionization of Na I against 
recombinations, of free; electrons and Na II, 

rn Na I = an Na II n e ' (4) 

and. assume, ail of the. sodium has a fixed- abundance relative to 
hydrogen- and. exists, only in neutral or singly ionized form, 

n Na = n Na. I + n Na II = {n e + n H I> <Na/H> (5) 


to arrive at the relation 


<Na/H> 


= n 


1 + (r/an e ) 

Na I + n 

n H I e- 


(6) 


The above equation may also be recognized for the relation of 
< Ca/H > with respect to n Ca jj. The accuracy of the ioniza- 
tion calculations, however, is severely limited by our incom- 
plete knowledge about conditions along any path through the 
interstellar medium, as well as some uncertainties in the 
ionization and recombination cross sections. Nevertheless it. 
is instructive to evaluate an estimate for total element abun- 
dances based on reasonable estimates for the relevant parame- 
ters, even if many of the choices seem to be rather arbitrary 
and hypothetical. 

A recent survey by Radhakrishnan, Murray, Lockhart and 
Whittle (19 71) of 21-cm emission and absorption, in the direc- 
tion of 35 extragalactic sources suggests that the hydrogen, is 
about equally divided between cool, dense clouds and a sur- 
rounding hotter medium. These observations, together with an 
interpretation by Hjellming, Gordon and Gordon (1969) : on con- 
ditions in the direction of three, pulsars;, suggest, that we 
could adopt the picture that, on the average, about 5% of any 
line of sight has roughly 5 H atoms, cm -9 at 6Q°K and the. re- 
maining 9 5% contains 0.26 atoms; cm -3 of hydrogen at 1000°K. 
Pressure balance is maintained in this model, and the overall 
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figure 8a . — Co rrelatlo n plot* ofi the Ny meas urements j \or those 
stars where we could obtain Ha I and Ca II column densities 
[of i their upper limits] by applying the doublet ratio meth- 
od to equivalent width measurement*, on, by Integrating the 
V line optical depths over velocity In the high resolution 
observations o l Hobbs {1969). The latter method otf obtain- 
ing Nncl I -bs considered more reliable than the former, and 
these better values are plotted with symbols having a bold 
border. The dashed lines correspond to our estimate {,or 
the average ratios o & Ha I and Ca II to H I o & 3.5 x io ~ 9 
and 1.5 x 10~ 9 , respectively. 


n H corresponds to 0.5 atoms cm -3 , which is a good representa- 
tive value for the more trustworthy measures of Na I, Ca II 
and H I column densities. For electrons produced from ioniza- 
tions of H I by X-rays or low energy cosmic rays , we expect to 
find 


n e an H° ’ S T° * 3 5 (7) 

if the ionization rate is constant everywhere (Habing 1969) . 
Densities of 0.042 and 0.026 electrons cm -3 for the cold and 
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F/gute. 86. — S&e F igute Sa capt-lon. 


hot portions, respectively, are consistent with equation (7) 
and the observed average of 0.03 electrons cm -3 for dispersion 
measures of pulsars at known distances (Davidson and Terzian 
1969, Lang 1971) . Our interpretation of the dispersion mea- 
sure integral includes an extra rms fluctuation term (Lerche 
1970) whose value corresponds to the expected variations of n e 
between the two regions. Habing (1969) estimates representa- 
tive values for r/a in the interstellar medium equal to 
2.1 (T/100 ) 0 * 7 for Na I and 0 .04 (T/100) 0 • 7 for Ca II. Substi- 
tuting these values into equation (6) we find < Na/H > = 2.3 

x 10" 7 and < Ca/H > = 6.8 x 10 " 9 is consistent with both the 

model and our averages for the observed < Na I/H I«> and 
< Ca II/H I > quoted earlier. These element abundance ratios 
are a factor of 7 and 230 below the cosmic abundance of 1.6 x 
10" 6 for both < Na/H > and < Ca/H > (Allen 1963). If we in- 
stead adopt a model for the interstellar medium closer to that 
suggested by Field, Goldsmith and Habing (1969) (3% of the vol- 

ume with nn = 30 atoms cm -3 and T = 60 °K surrounded by nn = 

0.1 atom cm -3 at 10 4o K) and apply the same analysis as before, 
we obtain the even lower values of <Na/H> = 1.0 x 10" 7 and 
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< Ca/H > = 4.0 x 10 " 9 . 

Figures 8a and 8b show the observed N Na i/ n H and N^ a n/1% 
ratios are generally lower for the Scorpius stars ( □ ) than 
for most of the others. This relative deficiency has been 
noted earlier for these stars by Buscombe and Kennedy (1968) 
in their comparisons of the interstellar lines with color ex- 
cesses. We previously suggested (in § IVa) that large Njj 
values measured for the Scorpius region may result from a rel- 
atively dense gas complex near these particular stars. We can 
surmise from Figure 7b that much of this gas is probably rath- 
er cool since, as an exception to the rule, the La Nh measure- 
ments are greater than the indications from 21-cm emission 
(i.e,, T 21-cm > -' L thus T^peak Tk~50°K). If we combine equa- 

tions (6) ana (7) with the T 0,7 dependence for F/a, we would 
expect the extra gas near the Scorpius stars, which is probab- 
ly cooler and denser than average, to exhibit larger 

< Na I/H I > and < Ca II/H I> ratios than usual — just the 
opposite to what is observed. It is possible, on the other 
hand, that the proximity of the stars to the gas may raise the 
ionizing radiation flux significantly above what is usually 
found in the galaxy. 

c) Correlation of H I with Dust 

A number of previous investigations have dealt with the re- 
lationship between interstellar hydrogen and interstellar dust. 
For a review of these observations see Kerr (1968) . All these 
studies have utilized 21-cm emission or absorption data along 
with a number of different techniques to obtain line-of-sight 
extinction. However, our poor correlation of 21-cm and L a 
column densities exposes the weakness of the assumption that 
the 21-cm radio observations of the gas and the optical obser- 
vations toward stars refer to the same regions of space. With 
the OAO L a absorption data and optical extinction data, one 
samples precisely the same regions of space so it is now pos- 
sible to make a more reliable study of the H I to dust corre- 
lation. 

Figure 9 shows the relationship between the OAO hydrogen 
column density and the E(B-V) color excess for all the stars 
included in Table 1. Lower case error symbols refer to B2 
stars while capital letters refer to stars of spectral type 
earlier than B2. The subscript e on a symbol indicates an 
emission line star for which the computed E (B-V) color excess 
is probably affected by the emission. 

In Figure 9 one can see a strong correlation between the hy- 
drogen column density and color excess. The B2 stars do not 
lie significantly above the earlier type stars, which indicates 
that for B2 stars it is reasonable to assume that the stellar 
L a line is contained within the core of the interstellar line 



LYMAN ALPHA ABSORPTION 


275 



Elguse 9. — A c ompaslhon o ^ mea*used Ny ves* u* 8-1/ colon exce** 
^os cull ol the *tas& ob* esved. Symbol* followed by an e 
am claAbl^led a* eml*t>lon line *tas*', we expect thels 
E(B-V) value * to be somewhat unsellable. Lowes ca*e essos 
symbol* Identify the 82 *tas* . The da*hed line sepse*ent * 
ous weighted avesage ofi 5 x U) Z1 atom * cm ~ 2 mag _1 los the 
ob*esved H I to E[B-V) satlo* . 


for most of the objects considered here. We found the situa- 
tion to be quite different for stars of spectral type B2.5 or 
later. For these stars, it was apparent from the line 
strengths in unreddened objects that stellar La was very 
strong. Because of these considerations we omitted all stars 
of spectral types later than B2 from this paper. 

The points in Figure 9 show a fair amount of scatter. Some 
of this is definitely observational while some is definitely 
real. For example, several Be stars deviate significantly 
from the general correlation. In the figure the scatter is 
largest for stars with small color excesses. Much of this 
scatter can be understood as small errors (0 I ?02) in the color 
excess determinations. In addition, the 0A0 column density 
measurements are less accurate for stars with small column den- 
sities due to the greater effect uncertainties in line blending 
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have on the L„ profiles. 

The data plotted in Figure 9 can be used to estimate the 
mass density ratio of H I to dust in the interstellar medium. 
Previous estimates from the radio data are summarized by Kerr 
(1968) . A mean of the numbers from various determinations 
(Kerr's Table 4) is ^ P H j/Pd us t > ~ 100, although the determi- 
nations vary from a low of 8 to a high of 300. Using the OAO 
L a data for stars other than Be stars and weighting each deter- 
mination by the color excess we obtain, 

<( N H /E(B-V) ) = 5 x 10 21 atoms cm -2 mag -1 . (8) 

To relate H I and dust densities it is necessary to assume 
a grain model. The radio investigations mentioned above gen- 
erally assumed a model consisting of a grain radius Ag = 3 x 
10“ 5 cm, a grain density Pg = 1 gm cm -3 , and an extinction 
efficiency factor Q = 2.0. With this very simple model plus 
the assumption that the ratio of selective to total visual ex- 
tinction Ry = 3.0, one can find from the relations given by 
Lilley (1955) and equation (8) that the OAO data imply 
K. P H l/ p dust^~ 50.’ 

However, it is possible to obtain an improved estimate by 
using the best interstellar grain model presently available. 

The theoretical mixture of grains of graphite, silicon carbide 
and silicates of Gilra (1971) has been shown to reproduce rea- 
sonably well the basic observational data (including ultravio- 
let extinction) on interstellar grains. Gilra (1971) finds on 
the basis of his mixture of grains that a space density of 8 x 
10~ 27 gm cm -3 is needed to produce 1 mag kpc -1 of extinction 
at the V filter. Again assuming Ry = 3.0 we find from equa- 
tion (8) that this mixture implies < Pjj i/Pdust)'- 100. 

The general appearance of the distribution of points in 
Figure 9 suggests the existence of a relatively well-defined . 
edge parallel to and above the dashed line which corresponds 
to the mean Nh/E(B-V). On the other hand, there is a signifi- 
cant straggling of the points below the line. We could specu- 
late that the upper edge (~2.5 times the number given in Equa- 
tion 8) corresponds to an intrinsic total hydrogen to dust 
ratio, and that lower points in the diagram show evidence of 
varying degrees of depletion of atomic hydrogen by either 
molecule formation or ionization. One outstanding example of 
possible depletion is t, Oph which is deficient of atomic hy- 
drogen by about a factor of 3 from the average ratio we de- 
rived. The visual interstellar spectrum of £ Oph is rich 
with molecular lines (Herbig 1968) , and it is possible molecu- 
lar hydrogen may also be abundant. McGee et. al. (196 3) ob- 
served a deficiency of 21-cm emission near S Oph which coin- 
cides with the Stromgren sphere for this 09.5 V star, but rep- 
resentative models for the size and density of the ionized re- 
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gion suggested by McGee et al. and Herbig imply the depletion 
of the observed Njj by ionization is negligible. 

Carruthers (1970b) has observed a molecular hydrogen column 
density of 1.3 x 10 20 molecules cm -2 in the direction of £ Per. 
The fractional abundance by weight of molecular compared to 
atomic hydrogen implied by this number is about 0.13 if one 
uses the OAO hydrogen column density for E, Per of 20 x 10 20 
atoms cm -2 . With such a low fractional abundance one might 
expect the N H /E(B-Vj relation for E, Per to be normal, and it 
is. An H 2 column density in the direction of £ Oph of 5 x io 20 
molecules cm -2 would explain the factor of 3 deficiency of H I. 
The implied difference in the H 2 abundance between these two 
stars might be explained by the fact that the far ultraviolet 
extinction in the direction of £ Oph is much greater than for 
E, Per (Bless and Savage 1972), and it is the far ultraviolet 
radiation field between 912 A and 1108 A that likely provides 
the destruction mechanism for the H 2 molecule (see Hollenbach 
et ad. CL971) for a discussion on the formation and destruction 
of H2) . 


V. CONCLUDING REMARKS 

The primary data of our study are the hydrogen column densi- 
ties toward the 69 stars listed in Table 1. The abstract sum- 
marizes a number of conclusions derived from a comparison of 
these column densities and other observables. We will not re- 
peat these conclusions here but will indicate some additional 
projects that could be undertaken. 

The comparison between 21-cm and L a column densities dis- 
cussed in § IVa would be more meaningful if one considered only 
the 21-cm observations with a small beam centered on each star 
and attempted to extract velocity components not appearing in 
the visual interstellar line data. For this project more star- 
centered radio observations are desirable. 

In estimating the total sodium and calcium abundance ratios, 
it should soon become possible to reduce some of the uncertain- 
ties. The sky background ultraviolet flux used by Habing 
(1969) to estimate F can hopefully be replaced by OAO-2 broad- 
band photometry observations of this radiation field. Also, 
according to Hobbs (private communication) , high rasolution ob- 
servations of the K line of Ca II should soon be available and 
will provide a more reliable Ca abundance estimate. 

Column 19 of Table 1 contains data from Duke (1951) , 

Stoeckly and Dressier (1964) and Wampler (1966) on the diffuse 
feature at 4430 A. It should now be possible to investigate 
if this feature is more closely correlated with the interstel- 
lar gas or the interstellar dust. C. C. Wu (private communica- 
tion) has undertaken an investigation of diffuse and semi-dif- 
fuse features in stars observed by the OAO. He has obtained 
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new observational data on the semi-diffuse features at 5780 A 
and 5>7/9'6' A and will investigate the relationship between the 
central, depths of 4430, 57 80 and 579 6 A versus the L a Njj, 
E(B-V) and E(2200 A-V) . In view of the work of Wu, we did not 
investigate the various correlations of diffuse features in 
this paper. 

There remain a number of additional stars B2 or earlier for 
which L a measurements can be made by OAO-2. We will try to 
observe these stars and fill some of the vacant areas in our 
sky coverage. Finally, as more high resolution prototype 
spectra become available, it should become possible to derive 
more accurate column densities for those OAO stars having 
large uncertainties. A rocket recently launched by the Kitt 
Peak National Observatory for the Princeton group obtained 
spectra for x Sco (BO V), k Sco (B1.5 III) and X Sco (B1.5 IV) 
Not only are these observations of excellent quality, but the 
actual column densities for the stars are relatively small so 
that reconstructions for small trial column densities will be 
more reliable. 

Support for the preparation of this paper was provided by 
contracts NAS 5-1348 (for BDS) and NSr-31-001-901 (for EBJ) 
from the U. S. National Aeronautics and Space Administration. 
We wish to thank Dr. A. D. Code for making the OAO data avail- 
able to us, and we are indebted to Drs . T. Houck, R. C. Bless, 
C. F. Lillie, A. Holm and M. Molnar for obtaining the observa- 
tional material. Important contributions were made by C. C. 
Wu, who helped reduce the raw OAO data to derive useful La 
profiles, and J. Pritchard, who computed Na I column densities 
from the observations of Hobbs (1969) . We are thankful to 
L. M,. Hobbs for his helpful comments and advice about his high 
resolution measurements. Dr. G. R. Carruthers kindly provided 
us with an original tracing of his y Ori spectrum. The Prince 
ton rocket spectra were obtained on NASA flights 4.176UG, 
4.226UG, and 4.271UG and were digitized by the microdensitome- 
ter at the Sacramento Peak Observatory. Finally, we wish to 
thank Mr. T. Matilsky for his help in deriving the intensity 
profiles of 6 and it Sco and £ Oph. 
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